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SOME APPLICATIONS OF A PROGRAMMING LANGUAGE 


I. SUMMARY 

APL is a programming language specifically designed for applied mathe- 
matics. It differs from other programming languages which are commonly used 
in that it is more concise, precise, and economical of symbols. In this paper 
is shown the APL language for some commonly used transcendental functions. 
Several improvements of APL are required for such things as program struc- 
turing, access to data, presentation of data, i.e., input/output operations and 
data operations. However, it is felt that this paper gives some indication of the 
usefulness of APL if one approaches it openmindedly. 

A Programming Language, APL is used for certain applications such as: 

1. Trigonometric functions and Solution of Kepler's Equation 

2. Harmonic Analysis 

3. Matrix Operations 

The areas of application are chosen primarily for their intrinsic interest 
to the Mission Trajectory Determination Branch and to indicate the usefulness 
of the language for Mission Trajectory Determination applications. 


II. TRIGONOMETRIC FUNCTIONS 

The APL language for computing the Arc cosine, Arc sine, Arc tangent 
(principal values), Arc tangent (quadrant oriented) Sine and Cosine are presented 
in Part A. The algorithm for all these computations except for the Arc tangent 
(principal values) are from reference 2. The algorithm for the Arc tangent 
(principal values) is presented in the Appendix. Part B is a program which 
computes the sine and cosine of several groups of data which are typical of some 
data used in the Mission Trajectory Determination Branch. The program along 
with the flow chart and the output of the program is presented. 
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PART A 


1. Arc Cosine, Arc Sine 

2. Arc tangent (principal values) 

3. Arc tangent (quadrant oriented) 

4. Sine, Cosine 

5. Solution of Kepler’s Equation 
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[37] Q*-Qx 57.29577951 

[38] ->-0 



Q = R Arcsincos A, 1 

Input: R, two element vector, if first element = 1 compute arcsin, if 

first element = 0 compute arcos; if second element = 1 answer 
in degrees, if second element = 0 answer in radians 

A, given sine or cosine 

Output: Q, computed angle in degrees or radians 
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B * ATAN A 

Input: A - tangent of unknown angle 

Output: B - value of angle in radians 



NO 


GET NEAREST STORED 
ANGLE WHOSE ARGUMENT 
IS SMALLER IN MAGNITUDE 
THAN PRESENT ANGLE 


t 

GET ARGUMENT OF 
ANGLE SUBTENDED 
BETWEEN TABLE ANGLE 
AND PRESENT ANGLE 



PUT GENERATED 
VALUE IN 

POLYNOMIAL SERIES 


B = ATAN A = 
SUM OF TABLE 
AND POLYNOMIAL 
VALUES 
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[26] J -*-J + 1 

[ 2 7 ] ->8 

[ 28 ] +0 


C •* A Qtan B, 3 

Computes angle C (quadrant oriented) given A = sin C and B = cos C 
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[33] +10 


R Sincos A, 4 


R = 2 element vector, A - angle, a , 0<a<277 

1st element = 1 compute sine, 1st element = 0 compute cosine, 2nd 
element = 1 , a in degrees, 2nd element = 0, a in radians 
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Solution of Kepler's Equation, 5 


VECCENTPICA NOML □ 3 V 
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Input is EC = eccentricity 

M = Mean anomaly 

Output is E = Eccentric anomaly 
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Input is EC = eccentricity 

E = Eccentric anomaly 

Output is F = True anomaly 
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PART B 


A program to compute the sine and cosine of some angles using APL 
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A 


B NTRIGROUT A 


Input: B is a vector whose first two components are the satellite's 

International Number, third component is the number of bulletins, fourth com- 
ponent is number of angles, and component 5 thru n are bulletin 
numbers . 

A is a matrix of angles . 

Output: Z is a matrix, whose first column is bulletin numbers, second 
column is angle Ai, third and fourth columns are sin Ai and cos Ai respectively, 
fifth column is Ai+1, and sixth and seventh columns are sin Ai+1 and cos Ai+1 
respectively. 
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C 3 5 ] * IMPLEMENTATION OF X. INCORRECT' 
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51 

*27.4651 

*0.4612082363 

0 . 8B72919297 

291 . 5981 

"0.9297886978 

0.3680937231 

52 

"46. 5842 

*0.7263851717 

0.6872878488 

229.388 

"0.7591349926 

*0.6509332279 

53 

"62.2016 

*0 .8845939995 

0.4663619427 

214.1932 

“0.561985219 

*0.8271472764 

54 

"81 . 0001 

*0.9876686147 

0.1 564327416 

195.3669 

"0 .2649991129 

”0.9642486631 

55 

263.9557 

*0 . 9944407755 

*0.1052973774 

179.5463 

0.007918476081 

*0.9999686511 

56 

249,864 

*0.9388781448 

*0.3442496798 

165.4586 

0,2510794908 

*0.9679664701 

57 

229.7032 

*0,7627044517 

*0.6467471867 

145.1242 

0. 5717994201 

*0.8203934583 

58 

217.4389 

0.6079150586 

"0. 7940020688 

132.6771 

0.7351855842 

*0.6778658874 

59 

183.4432 

*0.06005901175 

*0,9981948236 

99. 0378 

0. 9875849216 

*0.1570860433 

61 

146.6025 

0.5504443176 

“0 .8348718803 

60.6895 

0.8719795736 

0.4895422635 

62 

135.8983 

0.6959341056 

*0.7181056505 

49.8371 

0 . 764213813 

0.6449629851 

63 

127.7754 

0.7904180907 

"0,6125677477 

41 . 5232 

0.6629232615 

0 .7486873538 

64 

122 .5641 

0. 8427898097 

*0. 53B2428263 

36.6103 

0.5963691916 

0,8027102782 


BLTN 

>115 

SIN A 1 5 

COS A 15 

413 

SIN A 1 3 

COS A 1 3 

42 

271.1624 

"0.9997942111 

0.02028631552 

107.2344 

0.9551006558 

"0.2962815414 

44 

205.853 

*0.4360637343 

*0.8999157887 

43.5822 

0.6893945353 

0.7243860707 

45 

183. 7828 

*0.06597436079 

“0.9978213138 

19.6732 

0.336654853 

0.9416281236 

46 

147.2048 

0.5416377944 

"0.8406119813 

*16.5876 

*0.2854809615 

0.9583843879 

47 

106.8666 

0.9569727619 

*0.290177786 

"57.2938 

*0.8414523161 

0.5403313827 

48 

51.7166 

0,7849559019 

0,619551641 

113,2246 

*0.918966119 

"0.3943365086 

49 

0.151 

0.00263544412 

0.9999965316 

*165.011 

“0.2586335973 

*0.965975505 

50 

"28.2208 

*0.4728706768 

0 .881131845 

526.1112 

0.2400382962 

*0.970763426 

51 

"54.9302 

*0.8184526815 

0,5745739391 

583.1962 

*0.6844987525 

"0 . 7290140344 

52 

*93 . 1684 

*0 .9984713947 

*0 .05527083245 

458.776 

0.9882923796 

*0.1525718651 

53 

*124.4032 

*0 .8250819418 

"0.5650130906 

428.3864 

0.9296890811 

0.3683452417 

54 

*162.002 

*0.3089837982 

0.9510673092 

390.7338 

0.5110500707 

0.8595509411 

55 

527.9114 

0.2094240223 

*0.9778249275 

359.0926 

*0.01583645565 

0 .9998745961 

56 

499 . 728 

0. 646417003S 

*0.7629843129 

330.9172 

*0.48607306 

0.8739181798 

57 

459.4064 

0.9865539151 

0.1634361529 

290.2484 

“0.9382010068 

0 . 3460908605 

58 

434.8778 

0.9653716267 

0.2608785756 

265.3542 

"0.9967144476 

“0.08099568466 

59 

366.8864 

0 .11990116 

0.9927858281 

198.0756 

"0,3102716167 

*0.950647957 

61 

293.205 

*0 . 9191009617 

0.3940221214 

121.379 

0.853741699 

"0.5206967593 

62 

271.7966 

*0.9995084202 

0.03135144709 

99.6742 

0.9057792411 

0.1680455056 

63 

255.5508 

*0.9683692601 

0 .249521521 

83.0464 

0.9926445179 

0. 1210655C73 

64 

245.1282 

■ *0.9072511336 

*0.420589335 

73.2206 

0.9574233605 

0.2886875875 


BLTN 

414 

SIN A 14 

COSA 14 

>116 

SIN A 16 

COSA 1 6 

42 

190.1984 

*0 .1770572568 

*0.9842005555 

*7,8731 

"0.1369794933 

0.9905738828 

44 

124.7176 

0.8219691307 

*0.5695320461 

*73.0116 

"0.9563639364 

0. 2921780889 

45 

101.728 

0.9791235977 

"0.2032658104 

*94.1624 

"0.9973623199 

"0 . 07258369999 

46 

65.3086 

0.9085708912 

0.4177307076 

229.101 

*0.7558648963 

"0.6547276249 

47 

24.7874 

0 .4192524457 

0.9078697016 

188.9498 

*0.1555690371 

*0.9878250233 

48 

”30,754 

*0.5113530878 

0.8593707137 

134 . 1872 

0.7170663392 

"0.6970049281 

49 

*82.43 

*0 .9912846526 

0,1317373731 

82.732 

0.9919652526 

0.1265106106 

50 

248.9452 

*0.9332372472 

*0. 3592607026 

54.6132 

0.8152612294 

0.5790933699 

51 

222.133 

"0.670853858 

"0 .7415895803 

331.9934 

*0.4695732727 

0.8628935083 

52 

183.2038 

"0.05588772397 

*0.9984370552 

"9.5406 

"0.1657464504 

0 . 9861684024 

53 

151.9916 

0.4696010094 

"0.8828787558 

*40.796 

*0.6533941829 

0.7570178637 

54 

114.3668 

0.9109228839 

"0.4125766694 

*78.3672 

"0.9794599329 

0.2016386644 

55 

83.502 

0. 993575804 

0 ,1131685316 

252.3208 

*0.952771798 

"0.3036871997 

56 

55.3226 

0.8223685247 

0.568955194 

224.1334 

"0.6963313049 

"0.7177205017 

57 

14.8274 

0.2559080845 

0.9667011255 

183.9854 

"0.06950227415 

"0.9975817884 

58 

350.116 

*0 .1716539996 

0.9851573017 

1 59 .6396 

0.3479241642 

*0.937522687 

59 

282.481 

*0.9763677329 

0.2161158509 

91.2918 

0.9997458452 

"0.02254425317 

61 

207.292 

“0.4585254803 

*0.6886812645 

19.118 

0 . 3275147498 

0 .9448460734 

62 

185.7354 

*0.09993452248 

”0.9949940121 

"2.1422 

“0 .03737973299 

0.9993011304 

63 

169.2986 

0.1 856906257 

*0, 9B2608263 

*18.197 

*0.3122851799 

0.9499884107 

64 

159.1744 

0 .3555246141 

*0.9346669252 

*28.918 

“0.4835573996 

0.8753126562 


[ 
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IE. HARMONIC ANALYSIS 

The algorithm used here is presented in reference 3. The number of 
ordinates, i.e. , the number of divisions in the interval, can be either 12 or 24. 
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A 

V A-COMP HASH I 
Cl] +(COMPz 12)/ 5 

C 2 3 •*■( COMP= 24 )/4 3 

[33 'INCORRECT ORDINATE INPUT ' 

[4] ->0 

[5] 4-7p0 

[63 S-5p0 

[73 J-Hl 2 3 4 5 6] 

[83 K-Y[7 12 11 10 9 8} 

[9] 0-J+* 

[ 10 ] 

[111 C-U [ 1 2 3] 

[12] 0-i/[4 6 53 

[13] R-C+D 

[ 14 ] S-C-D 

[15] E-VZ 2 33 

[163 F«-V[6 S3 

[17] P-F + B 

[18] Q-E-F 

[19] L+-tf[2]+tf[3] 

[20] M^i?[ 2 ] -i?C 3 ] 

[213 0<?[1]+Q[2] 

[22] tf-«[l3-«[2] 

[233 4Cl3-(l*12)*(Stl]+L> 

[24] 4[2]«-(i*6)*( K[ l]t(0.5»( 3*0. 5)«5[2])*(0.S«5[ 3])) 

[25] 4[3Xl*6)*(S[l] + 0.5xtf) 

[26] 4[4]-(l*6)x(V'[i]-S[3]) 

[27] 4[5]«-(l*6)*(i?[l]-0.5xi) 

[28] 4[6]^<l*6)x(6r[l] + (*0.5x(3*O.5)x5C2]) + O.S*5[3]) 

[29] 4[7]-(l+12)x(5[i]-W) 

[30] *>(( f mi]-(+//U))>5S~6)/39 

[31] B[l]-(l* 6)x(P[4]*(0.5*P[l])+<0.S*(3*0.5)xp[2]>) 

[32] B [ 2]-*-(l + 12 ) x( 3*0. 5)*G 

[33] B[3]-(l*6)*(Ptl]-K[4]> 

[34] B[4]*-(l*12)«(3*0.5)*ff 

[35] B[ 5 ]—( 1 *6 ) * ( F[4] + (0.5xPCl3)+('"0. 5 *{ 3*0 . 5 ) *P[ 2 ] ) ) 

[36] 4-4 . B 

[37] -(< Kn2]-(B[l]+B[5] + (2xBC3]) + (3*0.5)x(F[2] + B[4]))))/4l 

[38] -0 

[39] ' CALCULATED A S INCORRECT ' 

[40] -0 

[41] * CALCULATED B 5 INCORRECT' 

[ 423 -*0 

[43] 4 — 1 3p 0 

[44] B-llpO 

[453 J-YZl 2 3 4 5 6 7 0 9 10 11 12] 

[46] *-Yt 13 24 23 22 21 20 19 18 1? 16 IS 14] 

[47] U-J + K 

[48] V-J-K 

[49] C\-U[l 2 3 4 S 6] 

[50] 01-B[7 12 11 10 9 8] 

[51] B-C1+D1 

[52] S-C1-D1 

[53] E\-*-V[ 2 3 4 5 6] 

[54] F 1-V[12 11 10 9 8] 

[55] P-Pl+Pl 
[563 Q+El-Fl 
[573 Bl-B[l 2 3] 

[58] J?2-B[4 6 5] 

[59] P-B1+B2 

[60] M-R1-R2 

[61] fll-Gtl 2] 

[ 623 C2-$[S 4] 

[633 C-Gl + <?2 

[64] H-Q1-Q2 

[65] E-l[2]i-Ll3] 

[663 F—L [ 2 ] -t[ 3 ] 

[67] <7-B[l]+B[2] 

[68] D-H [ 1 ]-ff[ 2 3 

[693 CO- 0 1 SINCOS 15 
[70] SI- 1 1 SINCOS 15 
[713 A[l]-(l*2*)*(l[l)+E) 

[72] 4[2]-(l*12>«< y[l]t(C0»5C2])+( ( 3 *0 . 5 ) *0 . 5 *S[ 3 ] ) + ( ( 1 ♦ 2 *0 . S ) *S[ 4 ] ) ♦ < 0 . 5 *SC 5] ) ♦ < SI*St 6 ] ) ) 

[73] 4[3]-(ltl2)x(S[l] + (0.Sx(3*0.S)*W[2])*(0.5x J w[3])) 

[74] 4[4]-( 1*12 )»(F[l] + (( 1*2*0. S)*(St2] + (-St4])+(-S[6]))> + <-S[S])) 

[75] 4[5] — (l*12)x(WC 13+0.5 *F ) 

[76] 4 [ 6 ]-<l* 12 )x( V[l]+(SI*B[2])+('0.5*( 3*0 . 5 ) *S[ 3 ] ) +( - ( 1 ♦ 2*0 . 5 ) *S[ 4 ] ) + { 0 . S»St 5 ] )+( COxSL 6 ] ) ) 

[77] 4[7]— {l*12)*(5[l } -Mt 3 ] > 

[783 4[B]-(l*12)x( FCl]+(-5J*SC2])+(“0.5«(3*0.5)«5C3])+(( 1 ♦ 2*0 . 5 )*S[ 4 ] ) +( 0 . 5xS£ S ] ) ♦< -0O*St 6 ] ) ) 

[79] 4[9]-(l*12)x(L[l]-(0.5xP)) 

[80] 4 [ 10]-( 1*12 )*( V[ 1 ]♦(-(! *2*0. 5 >*<S[2]+( -St «*]) + ( -Sf 6 ] ) ) ) ♦( -St 5 ] ) ) 

[81] 4[ll]-(l*12)*(Stl]+(~0.5x(3+0.5)xW£2])+(0.5xtf[3])) 

[82] 4[12]-( 1*12 )*< ni3 + ( -C0*5t2])+(0.5x( 3*0 . 5 ) *5[ 3 ] ) + ( -< 1 * 2*0 . 5 )*B[ 4 ])♦( 0 . 5 «S[ 5] )♦( -SI*Sl 6 3 ) ) 

[83] <4 [ 13 ]-( l*24)x(/¥[ 1 ] + (- F) ) 

CS»l] -(( !<Y[l3-( +//!>) )>5f'6>/39 

[853 B[ 1 ]— ( 1412 )*((SI x P[l3 )+( 0.5*P[2])*( ( 1*2*0. 5)*P[3])*( 0 . S *( 3 *0 . 5 ) xp[ 4 ] ) ♦( CO*Pl 5 ] ) + VZ 7 ] ) 

[86] fl[2]-(l*12>*<(0.5*Gtl])+< 0.5*( 3 * 0 . 5 ) *Gt 2 ] ) +<Jt 3 ] ) 

[87] fl[33-(l*12)*(F[2] + (-K[7])+(<l*2*0.5 WPti]*P[3] + ( -Pt5])))) 

[88] B[4]— (li24)x(3*0,5)*C 

[89] PC 5]-( 1*12 )x((CPxPCl] ) + ( 0.5xPC2]) + (-(l + 2*0.5 )*Pt 3] )+(’C.5*( 3 + 0.5 )*PC4]) + (SJ*Pr S])+y[73) 

[90] B[6]-(l*12>x(G[l] + (-<?U])) , „ , 

[91] S[7]*-(l*12)x((C0xPtl]) + CO.5xP[2])+(-(l*2 + O.5)xp[33)+<O.5x( 3 + 0 . 5 ) *PC 4 ] ) +{ SI* PC 5 3 ) + ( - V[ 7 ] ) ) 

[ 9 2 3 St B]-<1*24)*{3*0. 5 ) *0 

[93] B[9]-(l*12)x(y[7] + (-Pt2])+{<l*2 + 0.5)x(p[l] + Pt3] + ( -PC 5] ) ) ) ) 

[94] B[10]+-( 1*12 )x(( 0. 5xCC l]) + ("0. 5*0*0. 5)«Ct 2 ]>♦<?[ 3]) 

[95 3 fl[ 11 ]-{ 1*12) x((SIxPCl J) + ( “0.5xPt2] )♦(( 1«2*0.5 )*Pt 3] ) + (“0. 5*( 3*0. 5) xPr4 3) + (00xp£ 5]) ♦( -Ft 7])) 

[97] i(?ifn23-{((2xSI)x(Stl] + 8[H])) + {St2]+BtlO])*({2*O.5)x(SC3]*SC9]))*((3*O.5)x(BC4]+St0])>*((2xCP)x(p[ 

53+B[7]))+(2xBt6]))))>5P"6)/41 

[98] -0 
V 
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IV. MATRIX OPERATIONS 


At the time of this writing the matrix operations for APL had not been 
implemented. These programs, 1. INV (Matrix Inverse) and 2. TRANS (Matrix 
Transpose), show applications with existing operations. The Matrix Inverse 
should be especially noted for its conciseness when compared with other 
languages . 
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RB * Inv RA 

Input: RA, square matrix 

Output: RB, inverse of RA found by calculations in place 
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frwvCDlv 

V RB+INV RA ; RK ; RS ; RP ; RI 

[1] ->( ( 2 = ppRA)* = /pRA )p4 

[2] » INCORRECT MATRIX DIMENSIONS ' 

[3] ->0 

[4] /?£«-[/ pi?/l 
L 5 j 

[6] RP+\RK 

[73 RA <-RA [ ; ( i RS ) , 1 ] 

[ 8 ] RAl il+RSl+RSal 

[9] RI+{ \.RAl\RK',l])\[ / | RA [ t RK ; 1 ] 

[10] RPll ,RI1+RPIRI ,1] 

[11 j RALi.RI-, \RSl+RAlRI ,1 ; \i?S] 

[12] -*-((5Z7 — 10)> | RA [l;l])p21 

[13] 7i71[l;]«-/M[l;]*/M[l;l] 

[14] RA+RA -( (~RSa 1 ) *RA [ ; 1 ] ) <> ,x/f4[l ; 3 

[15] RA+-RAll+RS\ i RS ; ( 1 + i RS ) ,1] 

[16] RP+RPL1+RS\ \RS1 

[17] RK+RK- 1 

[18] -*■( RK > 0 ) p 8 

[19] RB+RA [ ',RP\ \i?S] 

[ 20 ] -*-0 

[21] 'NO INVERSE • 

[22] ->0 

V 
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VTi?i4il/5CQ]7 
V Z+TRANS Ail iJ 
Cl] 

C 2 ] Z+I [2 l]pO 

C 3 ] J -*-l 

C 4 ] ZU;>4[;«7] 

C 5 ] -*• ( XC 2 J -*-J + 1 ) / 4 

7 


4-«-5 6p 1 3 0 
A 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 


TRANS A 


1 

7 

13 

19 

25 

2 

8 

14 

20 

26 

3 

9 

15 

21 

27 

4 

10 

16 

22 

28 

5 

11 

17 

23 

29 

6 

12 

18 

24 

30 
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V. APPENDIX 


On A Method of Computing Tan* 1 x for Double Precision 

The evaluation of Arctangent x, in this method, is accomplished by a com- 
bination of "table look-up" and evaluation of a polynomial series. This elimi- 
nates the need for a subroutine for x < 1 and another for x > 1. The subroutine 
has stored a table of arguments and their corresponding angles (see table 1) . 

The "table look-up" part consists of obtaining the nearest stored angle whose 
argument, x is smaller in magnitude than the given argument. The subroutine 
then evaluates the argument of the angle subtended between the nearest table 
angle and the desired angle. The resulting argument is sent through a "table 
look-up" until the generated argument becomes smaller than the smallest non- 
zero argument in the table. The final generated argument is then inserted in 
the polynomial series . The value of the arc tangent is the sum of the values 
obtained by the "table look-up" process and by the polynomial series. 

Example 


Compute 6 = Arctan x given tan 6 , x = 1.7320508 

1. The closest value to x in the stored table (see table 1) is 
x = 1.000 giving a Q x = .7853 . . . radians. 

tan a - tan B 

2. Now from the relationship tan (a -B) = — j - 

H 1 + tan a • tan 

1 732058 —1 

one has tan (8-0,) ■= ^^7320508 = - 26794919 
so x j = .267949 . . . 

3 . Repeat step 1 using x x to find 6 2 = • 24497 . . . 

4. Repeat step 2 to find x 2 = . 01684 . . . 

Now x 2 is smaller than any non-zero table entry so x 2 is inserted in 
the polynomial series to yield .01682 . . . = d 3 . 

5. 9 then is equal to d\ + d 2 + 0 3 = 1.04719 .... 
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Table 1 


Tan 6 , x 
1 

.725 


.5 


1/3 


.25 


.1 


0 


0 *, Arc tan x 
.78539816 . . . 

. 62730819227 . . . 
.463647609 . . . 
.321750554397 . . . 
.244978663 . . . 
.099668652491 . . . 
0 


♦Note that 16 figures will be required to obtain the desired accuracy. 
Polynomial Series 

x(c 0 -x 2 (c i ~x 2 (c 2 -x 2 (c 3 -x 2 (c 4 -x 2 (c s -c 6 x 2 )))))) 

with 


1 _1 1 _ 1 1 _ 1 
c o~ l» c i g > c 2 “ g > c 3 ~ 7 ’ ° 4 9 * c 5 ~ ll ’ + 13 

The polynomial is "nested" to avoid the effects of multiplying by a power 
of x greater than 2 in any one operation. 

This method has been found to be accurate in a subroutine using up to 20 
significant figures in increments of 4 significant figures and for a range of 
~ 25.2 radians, i.e. , the approximate equivalent of 4 complete cycles. 
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